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The effects arising from the inherent continuous-mode nature of photonic pulses were poorly 
understood but significantly influence the performance of quantum devices employing photonic pulse 
interaction in nonlinear media. Such effects include the entanglement between the continuous wave- 
vector modes due to pulse interaction as well as the consequence of a finite system bandwidth. We 
present the first analysis on these effects for interactions between single-photon pulses, demonstrating 
their impact on the performance of quantum phase gates based on such process. Our study clarifies 
a realistic picture of this type of quantum devices. 



I. INTRODUCTION 

Deterministic photon-photon phase gate is a key build- 
ing block to construct circuits for the scalable all- 
optical quantum information processing. Cross-phase 
modulation (XPM) between slow pulses in media under 
electromagnetically induced transparency (EIT) condi- 
tions P, P (or with similar properties) is the main route 
towards such a gate. Considerable theoretical develop- 
ments (see, e.g., [3|-ll6|) as well as experimental studies 
(see, e.g., |17h23| ) have been undertaken to explore the 
EIT-based XPM and alternative approaches aiming at 
realizing photonic two-qubit gates. 

In optics-based quantum computing, a conditional 
phase of tt radians needs to be implemented by photon- 
photon phase gates. For achievin g su ch a large phase 
there are the proposals d, |^ 0, flolTll| of making pulses 
co-propagate, so that the generic weak pulse interaction 
could be compensated by the prolonged interaction time. 
Another requirement for an ideal phase gate is the uni- 
formness of conditional phase, as given by the mapping 
|l)i|l)2 — >■ e'^|l)i|l)2, i.e., the same phase 6 is induced 
for each mode fc of a continuous- mode photon in the state 
|1) = J^oo^^ C(^)^^(fc)|0), where ^{k) is the pulse pro- 
file in wave- vector space. The currently dominant un- 
derstanding is that a homogeneous phase could be pos- 
sible if the interaction between two pulses is averaged 
out by letting one pulse completely go through the other 

UMMMM- 

Since a photon is not a point particle, photon-photon 
interactions in nonlinear media should be modeled as in- 
teracting quantum fields of continuous modes. Using this 
picture, we show that the above-mentioned notions are 
generally invalid. The essential effects in realistic single- 
photon XPM are clarified here for the first time. 



II. CONTINUOUS-MODE EFFECTS IN 
PHOTON-PHOTON INTERACTIONS 

We first provide a theoretical framework for the XPM 
between photons. The interaction between two pulses in 
a medium of any type of atomic structure realizing EIT 
can be translated into that between two dark-state po- 
lariton fields ^/(z,i) = a/(A:)e*'=^dfc {I = 1,2) 



with [ai(fc),a](fc')] = 5ij6{k~k') [24|. By neglecting the 
pulse loss and deformation, as well as the possible self- 
phase modulation term which has no effect on photon- 
photon interactions [Tsl - fTsj . one has the following equa- 
tions of motion for the slowly varying and transversely 
well confined polariton fields 9, 13-15]: 

{dt + Vid,) i>i{z, t) = -iai{z, t)i>i{z, t) (1) 

where Vi are the pulse group velocities. The term 
ai{z,t) = J dz'A{z-z')'^l_^{z',t)'^3_i{z',t) could come 
from a general interaction potential A(z — z'). The pulse 
interaction in the experimentally studied XPM thus far 
J17-23], for instance, can be modeled by a contact po- 
tential A(z — z') = x^iz — z'), where x approximated by 
a real quantity is the nonlinear rate determined by the 
specific system parameters. The potential A(z — 2') con- 
sidered here acts instantaneously; see [l^l for a study on 
the non-instantaneous effects. From the field-theoretic 
viewpoint, Eq. ^ is obtained by the equation of motion 
ihdt-^^ = 6H/S^ for non-relativistic fields, where the 
Hamiltonian H = K + V consists of the kinetic term 
K = Y^f^^ J dzvi'^j{z)jW z'ifi{z) and the interaction 

term V = hjdzj dz'i'l{z)iil{z')A{z - z')*2(z')*i(^)- 
The pulse interaction would evolve the input state to 
t/(t)|l)i|l)2 = J dkj dk'C{k,k',t)d^k)b''{k')\0), where 
U{t) = Tcxp{-i J*dt'H{t')} (T denotes the time- 
ordering operation and h = 1 is adopted here- 
after). It is convenient to use the two-particle 
function, ij;{zi,Z2,t) = (0|\E'i(zi, t)*2(22, Ol^in) = 
^ J dk /d/fc'C(fc,fc',t)e*'=^ie^'='^^ to study the evolution 
of the initial state I'&in) = |l)i|l)2- To indicate how close 
a realistic XPM is to the ideal one |l)i|l)2 e^''|l)i|l)2 
for phase gates, we use the fidelity F . As figures of merit 
to characterize a realistic XPM, the fidelity F and con- 
ditional phase 9 are determined by the overlap [l^, ITsj 

^/Fe'^ = (*o|$out) = JdziJdz2%{zi,Z2,t)il;{zi,Z2,t), 

(2) 

where ipg [zi , Z2,t) is the two-particle function from the 
field operator equations {dt + Vidz)^i{z,t) = 0, corre- 
sponding to the freely evolved state |$o) = e~*^*]J&in). 
A similar formula in the discrete form is given in [26| . 
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The output state of a reahstic XPM, |$out) = 
J dk /dfc'C(fc, A;',t)at(fc)6'l'(fc')|0), is generally entangled 
between the wave-vector modes k and k' of the individ- 
ual photons. Such field modes entanglement was widely 
neglected in the previous researches on photon-photon 
gates, though its effect on XPM was conjectured in 
It is conceivable that the entanglement would lower the 
fidelity F. Yet, for clarifying the issue, a relation between 
the amount of entanglement generated in pulse interac- 
tion and the corresponding gate operation fidelity has 
to be found. Here we quantify the field mode entangle- 
ment with the linear entropy 5*^ = 1 — Trp| (pi are the 
reduced density matrices of the bipartite state l^^out)), 
which takes the following closed form: 

SL{t) = 1- J dzidz2dz3dzi^^{zi,Z2,t)'ip*{z3,Z2,t) 

X ■0(23,24,i)'0*(2l,^4,i)|; (3) 

see Appendix for the proof. Compared with the Schmidt 
decomposition method to characterize bi-photon entan- 
glement [27I [28j , this formula provides an exact measure 
of such entanglement. 

In reality the medium for pulse interaction carries 
a finite bandwidth Auj^, which is connected to the 
width of EIT transparency window [1, 0]- Hence, the 
equal-time commutator for the field operator 'i'i{z,t) ~ 
ai{k)e^''^dk after being imposed a cut-off by the 
bound of the wave- vector mode kg = Au}s/{2c) becomes 

{^.^{zi),^{z2)] = S^J{ks/n)amc{ksizl~~ Z2)) 

= (5,yC(zi-Z2), (4) 

where sinc(a;) = sin(x)/a;. Substituting the for- 
mal solution of ([1]), ^i{z,t) = exp{— i J*dt'ai(^z — 
Vi{t — — Vit,0), where ai{z) originates from 

the contact potential x^i^ ~ ^')^ vnto ip{zi, Z2,t) = 
{0\'i/i{zi,t)^2{z2,t)\^in) , while considering the commu- 
tator in (lU, yields the general two-particle function 



^izi,Z2,t) = (0|*i(zi)*2(4)l'fin> + (0|^ 



X / dtn--- / dtiC{Vrtn-l — Vrtn) ■ ■ ■ C{Vrti ~ Vrt2) 

Jo Jo 

XC{Z2 - z[- i;^ii)*i(z2 - Wri«)*2(4)l*in) (5) 



at the time t, where 



Vit and Vr = vi 



V2- 



The second term on the right side of Eq. ([S]) arises 
from photon-photon interaction. The norm of the two- 
particle function i/>(zi,Z2,i) is not preserved as a conse- 
quence of the deviation of the field operator commuta- 
tor C{z — z') from the delta function 5{z — z'), so the 
output two-photon function iplzi, Z2,t) has to be nor- 
malized prior to calculating the fidelity and linear en- 
tropy using Eqs. ^ and ([3]). The system bandwidth 
Aujs effectively results in a non-unitary evolution, though 



the wave- vector modes fc, k' of two photons are still 
continuous after imposing the cut-off k^. This could 
be understood by the restriction —kg < k,k' < kg 
on the matrix elements (fc, fc'|i/(t)|fc, /c') of the evolu- 
tion operator C/(t), causing the loss of the orthonor- 
mal relation for the matrix elements {k, k'\U{t)\k, k'); c.f. 
the inequality /^^^ dl j'^l^ dl'{k, k'\U\l, l'\U^q, q') 7^ 
dl JZo k'm, l'\U^\<l, q') = S{k - q)S{k' - 

q')- 



III. IMPACT ON PHOTONIC PHASE GATE 
PERFORMANCE 

An important situation we analyze here is the XPM 
between two co-propagating pulses [1, IH 0, [13, [HI ■ With 
the relative velocity Vr = 0, the two-particle function in 
Eq. ([5]) reduces to (see also Ref. [J) 

■lp{zi,Z2,t) = /i(zi - Vt)f2{z2 - Vt) + /i(z2 - vt) 

x/2(z2 - vt)smc{^{zi - Z2))(e** - 1), (6) 

where fi{z) = (0|5'j(z)|l), $ = x^^st/i^nc) and vi ^ 
V2 = V. Using the normalized form of this two-particle 
function, one will obtain from Eq. ([2]) the following rela- 
tions to determine the conditional phase and fidelity: 



tan 6* = 



Ci sin $ 



1 - Ci + Ci cos * ' 



F = 



l-4Ci(l-Ci)sin2f 
1 - 4(Ci - C2) sin^ f 



(7) 



(8) 



where Ci ^ J dZ^ J dZ2fl(Z^)h{Z2)\f2(Z2)?s\nc[ko{Zi- 

Z2)l C2 = /dZi/dZ2|/l(Z2)|2|/2(Z2)|2(sinc[A:o(Zi - 

■^2)])^- The dimensionless variables of the integrals are 
Zi — {zi — vt)/a, where a is the pulse size in medium; the 
system parameter is defined as fco = v/{2c)AuJs/Aujp, in 
proportion to the ratio of the system bandwidth Aujs 
to the pulse bandwidth Aujp (the reciprocal of the pulse 
duration). The factor v/(2c) in the parameter ko reflects 
the pulse compression in EIT media. 

In Fig. [1] obtained from Eq. ^ , one sees two differ- 
ent XPM patterns depending on the values of Ci and 
separated at the transitional point Ci = 0.5. This 
Ci value is the threshold above which the denominator 
1 — Ci + Ci cos on the right-hand side of Eq. ([7]) may 
become zero or negative with increasing <&, and it is uni- 
versal with respect to any pulse shape. As Ci decreases 
or increases across this point, the gate performance char- 
acterized by the achievable 6 values will undergo a tran- 
sition to the different pattern. For Ci < 0.5, the condi- 
tional phase 6 can never assume values greater than 7r/2 
no matter how large $ is. 

Without loss of generality, in our discussion below we 
consider the interaction between two identical pulses in 
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the Gaussian profile f{z) = {-^y-^)^ exp{— ^3-}. Then 

the boundary value Ci = 0.5 of the two patterns in Fig. 
1 corresponds to the system parameter fco ~ 2.5. As /cq 
increases (decreases) from this value, the coefficient Ci 
will become smaller (larger) into the respective pattern. 
Eqs. ([3]) and ([8]) yield the evolution of the field- mode 
entanglement and gate operation fidelity with the phase 
<I>, which is proportional to the pulse interaction time t. 
The results are illustrated in Fig. 2 for various ko values. 

The expected relation between and F manifests in 
the lower regime of Ci < 0.5 in Fig. 1, where a weaker 
entanglement, at any fixed $ value, is accompanied by a 
higher fidelity. In this regime the entanglement between 
the field modes totally vanishes in the limit fco — >■ oo. 
In the upper regime of Ci > 0.5, the field mode entan- 
glement disappears in the limit fcg — > 0. The fidelity 
value, however, goes down to zero at this point, for any 
non-zero $ value; c.f. Eq. ([5]) with a diverging C2 at 
ko — > 0. The non-unitary evolution due to a finite sys- 
tem bandwidth Aujs accounts for this phenomenon. Note 
that, as the consequence of non- unitary evolution, both 
factors /i, /2 in the second term of Eq. ^ carry the 
same variable Z2 — vt, thus deviating from the ideal out- 
put two-particle function fi{zi — vt)f2{z2 — vt)e^^ in the 
neighborhood of fco — 0. The XPM between two evenly 
distributed square pulses shows the same characteristics 
except for the different relations between the Ci and fco 
values, indicating that the effects described above differ, 
by nature, from the inhomogeneity in pulse interactions. 

An interesting observation from Fig. 2 is the feature 
that, around the transitional point fco ~ 2.5 of the two 
operation patterns, the plots of Sl, e.g., the thin solid 
line, reaches the highest value at some $ < tt. Increasing 
$ beyond the point leads to a continual distortion of the 
two-photon state without increasing its entanglement. Ir- 
respective of the pulse shapes, the linear entropy plateau 
at the transitional point is correlated to the fidelity valley 
like that in Fig. 3, which goes down to = at fco « 2.5 
and $ = TT. 

Furthermore, we comment on the notion that pulses 
moving with matched group velocities make large con- 
ditional phase possible In the previous theoretical 








FIG. 1: (color online) Relation between the conditional phase 
9 and the XPM phase $ for various Ci values from the lower 
to the upper: Ci = 0.20, 0.36, 0.45, 0.50, 0.55, 0.63, 0.78 and 
0.98. The line 6 = ^"1? corresponding to Ci = 0.5 separates 
the different operation patterns. 




FIG. 2: (color online) Linear entropy Sl (left) and fidelity 
F (right) plotted vs $ for various fco values: fco = 0.5 (thin 
dashed line), fco = 1.0 (thick solid line), fco — 2.5 (thin solid 
line), fco — 5.0 (thick short dashed line), and fco = 10.0 (thick 
long dashed line). The two pulses are in the identical Gaussian 
profile. 

studies, the variables zi and Z2 of the two-particle func- 
tion in ^ were often mixed up with the pulse-center 
coordinates, and then the phase $ would be regarded 
as the conditional phase 6 for the reason that the over- 
lapped pulse centers with zi = Z2 = z could lead to the 
ideal output /i(z — tii)/2(z — fi)e** from Eq. In fact, 
the variable AuJs{zi — Z2)/(2c) of the sine function in ^ 
assumes any value even if the two pulses co-propagate, 
because Zi are the field coordinates over the whole z axis 
rather than those of the pulse centers. The conditional 
phase 9 should be determined by Eq. ^ giving its rela- 
tion with the XPM phase $ in Fig. 1. The conditional 
phase value could reach tt in the upper regime of Fig. 1, 
where the gate operation fidelity is, however, rather low. 
A high fidelity is possible only in the lowest region in 
Fig. 1, where the state evolution is close to unitary but 
the peak of the conditional phase 9 is vanishing. The 
vanishing conditional phase in the regime of near uni- 
tary evolution also exists in case of interaction between 
co-propagating single photon and coherent state; see p^ . 

The trade-off between conditional phase and fidelity in 
XPM between co-propagating photons is also discussed 
recently in [26| . where a finite mode approximation for 
pulses is adopted. One problem with the finite mode ap- 
proximation is its underestimation of XPM intensity — 
compared with the contribution from the second term in 
our Eq. ([S]), the term arising from XPM, the discrete sum 




FIG. 3: (color online) Gate operation fidelity _F as a function 
of the parameters fco and The parameter ranges are 0.1 < 
fco < 8 and < "1> < tt. It is a 3D view for the right plot of 
Fig. 2. 
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F 

i.Or - - - > B a fixed value from the contact potential x'^i^ ~ ^')- 

i B This result is also true to the XPM between single pho- 

« ■ ■ ■ ■ coherent state [1^1 ■ The time-dependent two- 
particle functions for such idealized unitary evolution un- 

1 ^ ♦ der a general interaction potential A(z — z^ take the 

<I>=0.1 ^ form /i(zi - wit)/2(z2 - W2t)e'*(^i^^^'*) The 

* * * field mode entanglement exhibited by the possibly non- 
factorisable $(zi,Z2,t) with respect to its spatial vari- 

* ables is therefore the main factor that determines the 

^ gate operation fidelity in a regime of approximately uni- 

* * * tary evolution, which could be realized under the condi- 
0.0005 0.0007 0.0009 0.0011 0.0013 0.0015 ^ tiou > AWp. 



♦ <I>=0.2 

* 'I>=0.3 



FIG. 4: (color online) Fidelity evolution in head-on collision 
between two identical single-photon pulses of Gaussian pro- 
file. The two pulses are initially separated by a distance of 
I = 10a, where a is the pulse size in the medium, and run 
toward each other at a relative velocity Vr — 2\vi\ = lO^cr per 
unit time (here we adopt a scaled time unit). The pulses 
completely overlap at t = 0.001. The system parameter 
fco = \vi\/{2c)ALOs/Aojp is 0.001. The XPM phase is defined 
as <I> = xAcjsZ/(7rurc). 

with the corresponding term in [26] (the second term in 
Eq. (22) of [2^) contributes much less significantly to a 
coefficient similar to Ci in this paper, limiting its value 
to less than 0.5. Then only the lower regime in our Fig. 
1 can be obtained in the finite mode approach. For two 
ultraslow pulses the small ratio v/c in the parameter ko 
makes the upper regime in Fig. 1 more relevant. The 
effect of non-unitary evolution of quantum states dom- 
inating in this regime is beyond the description by the 
finite mode approximation. 

Next we examine the XPM between pulses colliding 
head-on d, H, [l^, [13) via a contact potential xS{z — z'). 
For two ultraslow pulses satisfying Vr/{2c)Aujst <C 1, 
there is the approximation C{vrtk-i — Vrtk) ~ Aoj^/ (27rc) 
in ([5|), so the infinite sum in Eq. ^ can be approxi- 
mated by a closed form in this regime. Substituting the 
normalized two-particle function into ([2]), we obtain the 
fidelity evolution in the course of pulse interaction; see 
the example in Fig. 4. It shows that the fidelity value 
will decline once a pulse touches the other and stabilize 
again after they pass through each other. The stronger 
the interaction (indicated by the $ values) is, the lower 
the fidelity will become after collision. Against the intu- 
itive notion that an averaged interaction on pulses could 
generate a uniform conditional phase 9, a realistic XPM 
between pulses of a non-zero relative velocity can be far 
away from the ideal process |l)i|l)2 e*^|l)i|l)2- 

An ideal performance of phase gates based on the 
XPM described above exists in the limit of an infinite 
system bandwidth Awg. In this limit the two-particle 
function after two pulses completely going through each 
other is fi{zi — vt)f2{z2 + vt)e^^ , where 9 — x/vr is 



IV. CONCLUSION 

We have illustrated the effects of continuous field mode 
entanglement arising from pulse interaction and non- 
unitary evolution caused by finite system bandwidth, 
which drastically impair photon-photon phase gate per- 
formance. These effects induce more complexity in XPM 
than what was previously understood. Due to their possi- 
ble existence in any device working with quantum objects 
of continuous degrees of freedom, the proper handling of 
the effects could be a major concern in quantum technol- 
ogy 



Appendix 

We provide a brief derivation for the linear en- 
tropy formula in Eq. The elements of one re- 
duced density matrix pi for a general bipartite state 
fdkj dk'Cik,k',t)a^k)b^k')\0) can be obta,ined by the 
following [231 (its discrete form is given in |30j): 

pi{k,k' ,t) ^ J dq({k,q,t)Cik' ,q,t) ^ ^ J dzi J dz2 

X J dz3{i^{zi,Z2,t)r{z3,Z2,t)e-'''''e''''''}. 

The matrix elements of Pi can thus be obtained 
by substituting the above into p\{k,k' ,t) = 
J dqpi{k,q,t)pi{q,k' ,t), which leads to a closed form of 
the linear entropy SL{t) = 1 — J dkp1(k,k,t). 
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